Bias in ligation-based small RNA sequencing library construction is determined by adaptor and RNA structure.

High-throughput sequencing (HTS) has become a powerful tool for the detection of and sequence characterization of microRNAs (miRNA) and other small RNAs (sRNA). Unfortunately, the use of HTS data to determine the relative quantity of different miRNAs in a sample has been shown to be inconsistent with quantitative PCR and Northern Blot results. Several recent studies have concluded that the major contributor to this inconsistency is bias introduced during the construction of sRNA libraries for HTS and that the bias is primarily derived from the adaptor ligation steps, specifically where single stranded adaptors are sequentially ligated to the 3' and 5'-end of sRNAs using T4 RNA ligases. In this study we investigated the effects of ligation bias by using a pool of randomized ligation substrates, defined mixtures of miRNA sequences and several combinations of adaptors in HTS library construction. We show that like the 3' adaptor ligation step, the 5' adaptor ligation is also biased, not because of primary sequence, but instead due to secondary structures of the two ligation substrates. We find that multiple secondary structural factors influence final representation in HTS results. Our results provide insight about the nature of ligation bias and allowed us to design adaptors that reduce ligation bias and produce HTS results that more accurately reflect the actual concentrations of miRNAs in the defined starting material.

RNAi triggered by specialized machinery silences developmental genes and retrotransposons.

RNA interference (RNAi) is a conserved mechanism in which small interfering RNAs (siRNAs) guide the degradation of cognate RNAs, but also promote heterochromatin assembly at repetitive DNA elements such as centromeric repeats. However, the full extent of RNAi functions and its endogenous targets have not been explored. Here we show that, in the fission yeast Schizosaccharomyces pombe, RNAi and heterochromatin factors cooperate to silence diverse loci, including sexual differentiation genes, genes encoding transmembrane proteins, and retrotransposons that are also targeted by the exosome RNA degradation machinery. In the absence of the exosome, transcripts are processed preferentially by the RNAi machinery, revealing siRNA clusters and a corresponding increase in heterochromatin modifications across large domains containing genes and retrotransposons. We show that the generation of siRNAs and heterochromatin assembly by RNAi is triggered by a mechanism involving the canonical poly(A) polymerase Pla1 and an associated RNA surveillance factor Red1, which also activate the exosome. Notably, siRNA production and heterochromatin modifications at these target loci are regulated by environmental growth conditions, and by developmental signals that induce gene expression during sexual differentiation. Our analyses uncover an interaction between RNAi and the exosome that is conserved in Drosophila, and show that differentiation signals modulate RNAi silencing to regulate developmental genes.

Small RNA expression profiling by high-throughput sequencing: implications of enzymatic manipulation.

Eukaryotic regulatory small RNAs (sRNAs) play significant roles in many fundamental cellular processes. As such, they have emerged as useful biomarkers for diseases and cell differentiation states. sRNA-based biomarkers outperform traditional messenger RNA-based biomarkers by testing fewer targets with greater accuracy and providing earlier detection for disease states. Therefore, expression profiling of sRNAs is fundamentally important to further advance the understanding of biological processes, as well as diagnosis and treatment of diseases. High-throughput sequencing (HTS) is a powerful approach for both sRNA discovery and expression profiling. Here, we discuss the general considerations for sRNA-based HTS profiling methods from RNA preparation to sequencing library construction, with a focus on the causes of systematic error. By examining the enzymatic manipulation steps of sRNA expression profiling, this paper aims to demystify current HTS-based sRNA profiling approaches and to aid researchers in the informed design and interpretation of profiling experiments.

Structural bias in T4 RNA ligase-mediated 3'-adapter ligation.

T4 RNA ligases are commonly used to attach adapters to RNAs, but large differences in ligation efficiency make detection and quantitation problematic. We developed a ligation selection strategy using random RNAs in combination with high-throughput sequencing to gain insight into the differences in efficiency of ligating pre-adenylated DNA adapters to RNA 3'-ends. After analyzing biases in RNA sequence, secondary structure and RNA-adapter cofold structure, we conclude that T4 RNA ligases do not show significant primary sequence preference in RNA substrates, but are biased against structural features within RNAs and adapters. Specifically, RNAs with less than three unstructured nucleotides at the 3'-end and RNAs that are predicted to cofold with an adapter in unfavorable structures are likely to be poorly ligated. The effect of RNA-adapter cofold structures on ligation is supported by experiments where the ligation efficiency of specific miRNAs was changed by designing adapters to alter cofold structure. In addition, we show that using adapters with randomized regions results in higher ligation efficiency and reduced ligation bias. We propose that using randomized adapters may improve RNA representation in experiments that include a 3'-adapter ligation step.

T4 RNA ligase 2 truncated active site mutants: improved tools for RNA analysis.

BACKGROUND: T4 RNA ligases 1 and 2 are useful tools for RNA analysis. Their use upstream of RNA analyses such as high-throughput RNA sequencing and microarrays has recently increased their importance. The truncated form of T4 RNA ligase 2, comprising amino acids 1-249 (T4 Rnl2tr), is an attractive tool for attachment of adapters or labels to RNA 3'-ends. Compared to T4 RNA ligase 1, T4 Rnl2tr has a decreased ability to ligate 5'-PO4 ends in single-stranded RNA ligations, and compared to the full-length T4 Rnl2, the T4 Rnl2tr has an increased activity for joining 5'-adenylated adapters to RNA 3'-ends. The combination of these properties allows adapter attachment to RNA 3'-ends with reduced circularization and concatemerization of substrate RNA. RESULTS: With the aim of further reducing unwanted side ligation products, we substituted active site residues, known to be important for adenylyltransferase steps of the ligation reaction, in the context of T4 Rnl2tr. We characterized the variant ligases for the formation of unwanted ligation side products and for activity in the strand-joining reaction. CONCLUSIONS: Our data demonstrate that lysine 227 is a key residue facilitating adenylyl transfer from adenylated ligation donor substrates to the ligase. This reversal of the second step of the ligation reaction correlates with the formation of unwanted ligation products. Thus, T4 Rn2tr mutants containing the K227Q mutation are useful for reducing undesired ligation products. We furthermore report optimal conditions for the use of these improved T4 Rnl2tr variants.

Linear group II intron RNAs can retrohome in eukaryotes and may use nonhomologous end-joining for cDNA ligation.

Mobile group II introns retrohome by an RNP-based mechanism in which the excised intron lariat RNA fully reverse splices into a DNA site via 2 sequential transesterification reactions and is reverse transcribed by the associated intron-encoded protein. However, linear group II intron RNAs, which can arise by either hydrolytic splicing or debranching of lariat RNA, cannot carry out both reverse-splicing steps and were thus expected to be immobile. Here, we used facile microinjection assays in 2 eukaryotic systems, Xenopus laevis oocyte nuclei and Drosophila melanogaster embryos, to show that group II intron RNPs containing linear intron RNA can retrohome by carrying out the first step of reverse splicing into a DNA site, thereby ligating the 3' end of the intron RNA to the 5' end of the downstream exon DNA. The attached linear intron RNA is then reverse transcribed, yielding an intron cDNA whose free end is linked to the upstream exon DNA. Some of these retrohoming events result in the precise insertion of full-length intron. Most, however, yield aberrant 5' junctions with 5' exon resections, 5' intron truncations, and/or extra nucleotide residues, hallmarks of nonhomologous end-joining. Our findings reveal a mobility mechanism for linear group II intron RNAs, show how group II introns can co-opt different DNA repair pathways for retrohoming, and suggest that linear group II intron RNAs might be used for site-specific DNA integration in gene targeting.

EcI5, a group IIB intron with high retrohoming frequency: DNA target site recognition and use in gene targeting.

We find that group II intron EcI5, a subclass CL/IIB1 intron from an Escherichia coli virulence plasmid, is highly active in retrohoming in E. coli. Both full-length EcI5 and an EcI5-DeltaORF intron with the intron-encoded protein expressed separately from the same donor plasmid retrohome into a recipient plasmid target site at substantially higher frequencies than do similarly configured Lactococcus lactis Ll.LtrB introns. A comprehensive view of DNA target site recognition by EcI5 was obtained from selection experiments with donor and recipient plasmid libraries in which different recognition elements were randomized. These experiments suggest that EcI5, like other mobile group II introns, recognizes DNA target sequences by using both the intron-encoded protein and base-pairing of the intron RNA, with the latter involving EBS1, EBS2, and EBS3 sequences characteristic of class IIB introns. The intron-encoded protein appears to recognize a small number of bases flanking those recognized by the intron RNA, but their identity is different than in previously characterized group II introns. A computer algorithm based on the empirically determined DNA recognition rules enabled retargeting of EcI5 to integrate specifically at 10 different sites in the chromosomal lacZ gene at frequencies up to 98% without selection. Our findings provide insight into modes of DNA target site recognition used by mobile group II introns. More generally, they show how the diversity of mobile group II introns can be exploited to provide a large variety of different target specificities and potentially other useful properties for gene targeting.

Group II intron-based gene targeting reactions in eukaryotes.

BACKGROUND: Mobile group II introns insert site-specifically into DNA target sites by a mechanism termed retrohoming in which the excised intron RNA reverse splices into a DNA strand and is reverse transcribed by the intron-encoded protein. Retrohoming is mediated by a ribonucleoprotein particle that contains the intron-encoded protein and excised intron RNA, with target specificity determined largely by base pairing of the intron RNA to the DNA target sequence. This feature enabled the development of mobile group II introns into bacterial gene targeting vectors ("targetrons") with programmable target specificity. Thus far, however, efficient group II intron-based gene targeting reactions have not been demonstrated in eukaryotes. METHODOLOGY/PRINCIPAL FINDINGS: By using a plasmid-based Xenopus laevis oocyte microinjection assay, we show that group II intron RNPs can integrate efficiently into target DNAs in a eukaryotic nucleus, but the reaction is limited by low Mg(2+) concentrations. By supplying additional Mg(2+), site-specific integration occurs in up to 38% of plasmid target sites. The integration products isolated from X. laevis nuclei are sensitive to restriction enzymes specific for double-stranded DNA, indicating second-strand synthesis via host enzymes. We also show that group II intron RNPs containing either lariat or linear intron RNA can introduce a double-strand break into a plasmid target site, thereby stimulating homologous recombination with a co-transformed DNA fragment at frequencies up to 4.8% of target sites. Chromatinization of the target DNA inhibits both types of targeting reactions, presumably by impeding RNP access. However, by using similar RNP microinjection methods, we show efficient Mg(2+)-dependent group II intron integration into plasmid target sites in zebrafish (Danio rerio) embryos and into plasmid and chromosomal target sites in Drosophila melanogster embryos, indicating that DNA replication can mitigate effects of chromatinization. CONCLUSIONS/SIGNIFICANCE: Our results provide an experimental foundation for the development of group II intron-based gene targeting methods for higher organisms.